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Abstract: We report here on two years of timing of 168 pulsars using the Parkes radio telescope. The 
vast majority of these pulsars have spin-down luminosities in excess of 10^** erg s"^ and are prime target 
candidates to be detected in gamma-rays by the Fermi Gamma-Ray Space Telescope. We provide the 
ephemerides for the ten pulsars being timed at Parkes which have been detected by Fermi in its first year 
of operation. These ephemerides, in conjunction with the publicly available photon list, can be used to 
generate gamma-ray profiles from the Fermi archive. We will make the ephemerides of any pulsars of 
interest available to the community upon request. In addition to the timing ephemerides, we present the 
parameters for 14 glitches which have occurred in 13 pulsars, seven of which have no previously known 
glitch history. The Parkes timing programme, in conjunction with Fermi observations, is expected to 
continue for at least the next four years. 
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1 Introduction 

The Fermi Gamma-Ray Space Telescope, launched in 
2008 June, has the study of pulsars at gamma-ray 
wavelengths as one of its key science projects. Fermi 
follows in the footsteps of its predecessor, the Comp- 
ton Gamma-Ray Observatory (CGRO) which was ac- 
tive during the 1990s. CGRO detected seven gamma- 
ray pulsars (including the radio-quiet pulsar Geminga) 
with three other possible detections (including a mil- 
hsecond pulsar). CGRO relied on ephemerides for 
the pulsars it observed in order to correctly phase tag 
the received gamma-ray photons to produce a light 
curve. Three of the pulsars detected by CGRO (PSRs 
B1706-44, B1055-52 and B1509-58) used data ob- 
tained with the Parkes radio telescope (Thompson et 
al. 1992, Fierro et al. 1993, Ulmer et al. 1993, John- 
ston et al. 1995) . 

The Large Area Telescope (LAT) on board Fermi 
(Atwood et al. 2009) has a field of view which is sub- 
stantially bigger than that of the Energetic Gamma- 
Ray Experiment Telescope (EGRET) aboard CGRO, 
has a sensitivity more than an order of magnitude 
greater, far superior energy resolution and a positional 
accuracy measured in arcmin rather than degrees. Fermi 
began sky-survey observations on 2008 August 11 and 
in survey mode observes the entire sky every three 
hours. The characteristics of the photons received by 
Fermi remained proprietary for the first year of oper- 



ation but full data release occurred in 2009 August. 

In principle, therefore, Fermi can produce light 
curves for every pulsar in the sky and update the pro- 
file on a daily basis. To produce these profiles for all 
but a handful of the brightest pulsars, generally re- 
quires an accurate pulsar ephemeris so that the in- 
coming photons can be correctly phase tagged. Prior 
to the launch of Fermi, the worldwide radio pulsar 
timing community formalised a comprehensive pulsar 
monitoring campaign using the Parkes, Lovell, Nangay, 
Green Bank and Arecibo telescopes (Smith et al. 2008) . 
The time pressure on these large telescopes means that 
it is not possible to time all of the ~2000 known pul- 
sars. The timing list created by Smith et al. (2008) 
contained 224 pulsars with a spin-down energy loss 
rate, E, larger than lO^** erg s"^. Of these 224 pul- 
sars, 156 are timed at Parkes. 

One of the major issues facing the radio timing 
campaign is that many of the pulsars potentially de- 
tectable by Fermi are young, high E objects. These 
pulsars are far from perfect clocks and suffer fr om a 
high degree of timing noise (|Hobbs et al ] |2004l 1. In 
order to track the timing noise to a level of accuracy 
sufficient to time-tag the photons to 10~^ — 10~^ of the 
rotational period, we need to observe the entire sample 
approximately once per month. 

Fermi has been a spectacular success since its launch. 
It has discovered 16 g amma-ray pulsars t hrough 'blind' 
periodicity searches (|Abdo et al.ll2009dh of which two 
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have subsequently been detected in the radio (Camilo 
et al. 2009). In addition it has produce d high qual- 
ity pr ofiles for 8 millisecond radio pulsars (jAbdo et alj 
l2009j ) and more than 20 'norma l' (non-millisecond ) 
pulsars, includin g the Crab pulsar ( Abdo et al.ll20093 ). 
the Vela pulsar (|Abdo et a l. 2009h ) ). the recently dis- 
covere d PSR J1028-5819 (|Keith et al.ll2008l:lAbdo et al.l 
l2009al l and 6 pulsars with moderate E (Weltevrede et 
al. 2009). All these results h ave been tabulated in the 
gamma-ray pulsar catalogue (jAbdo et al.ll20093 ). 

The radio data obtained through the pulsar tim- 
ing program has a number of other applications, in- 
cluding a statistical determination of the alignment 
of th e magnetic and rotational axis as a function of 
time ([Weltevrede fc J ohnston" 2008al 'l and examination 
of the p olarisation ch aracteristi cs of hi ghly energetic 
pulsars ([Weltevrede fc Joh nston 2008bY In addition, 
long term r adio timing shows tha t a number of pul- 
sars glitch llShemar fc Lvnel 1 19961 : IWang et all l2000l : 
Ijanssen fc Stapperd 120061 ) and the glitch parameters 
and time between events can be used to determine the 
interior structure of neutron stars (jRuderman et all 
Il998l l. 

Of the 'normal' radio pulsars detected by Fermi, 
Parkes is responsible for timing ten of them. In this 
paper we describe the Parkes timing campaign for the 
Fermi mission, including the observational and data 
analysis details. We provide the timing ephemerides 
for, and a brief discussion of, the 10 pulsars detected by 
Fermi. We tabulate the glitches seen in the monitored 
pulsars during the two years of the timing program. 
We briefly discuss the implications of the results and 
the next steps for the Fermi mission. 

2 Observations and Data Anal- 
ysis 

Observations in support of the Fermi mission com- 
menced in 2007 February. All observations were car- 
ried out using the 64-m radio telescope in Parkes, NSW, 
Australia. Generally, each observing session used an 
observing frequency near 1.4 GHz and lasted 24 hr 
with the observing sessions separated by approximately 
4 weeks. However, once every six months we had an 
extended observing session and data were obtained at 
both 3.1 and 0.7 GHz in order to monitor long-term 
dispersion measure variations. Each pulsar is typically 
observed for only a few minutes, sufficient to reach a 
signal-to-noise ratio greater than 5. A total of 156 
pulsars from the Smith et al. (2008) list are timed at 
Parkes; their basic parameters are given in Table [T] 
A further 12 pulsars are also timed; these are listed 
separately in Table [2] 

At 1.4 GHz, we used the centre beam of th e multi- 
beam receiver ijStavelev-Smith fc et al.l Il996l ) with a 
total bandwidth of 256 MHz. The noise-equivalent flux 
density of the system is ~35 Jy on cold sky. For the 
3.1 GHz and 0.7 GHz observing, data were reco rded si- 
mult aneously using the dual 10/50-cm receiver (|Granet fc 
120051 1 with a usable bandwidth of 1024 and 40 MHz re- 
spectively. The noise-equivalent flux density is ~49 Jy 
at 3.1 GHz and 70 Jy at 0.7 GHz. Many pulsars in our 



sample are located in the Galactic plane where the sky 
temperature can be up to -^50 K at 1.4 GHz and sev- 
eral hundred K at 0.7 GHz, substantially raising the 
noise- equivalent flux density on these sources. 

The downconverted signals from each polarisation 
channel of the linear feeds were fed in a digital fil- 
terbank, designed and built at the Australia Telescope 
National Facility. The hardware converts the analogue 
voltages into digital signals and produces a filterbank 
output consisting of 1024 frequency channels for each 
of 1024 phase bins across the pulse period for the two 
auto correlations and real and imaginary parts of the 
cross correlations of the feed probes. Data were folded 
at the topocentric pulse period and accumulated for 
30 s and then dumped to disk. A calibration signal, 
injected into the feed at an angle of 45° to the probes, 
was recorded on a regular basis. This signal is used 
to determine the relative gain of the two polarisation 
channels and the phase between them. 

Data analysis was carried out using the psrchive 
package (| Hot an et al.l |2004| ) . In brief, after removal 
of interference in both the frequency and time do- 
mains, the data are gain and polarisation calibrated 
and summed in frequency and time to produce a pulse 
profile for each pulsar observed. 

3 Pulsar Timing 

In order to produce a time-of-arrival (TOA) the resul- 
tant pulse profile is cross correlated with a high signal- 
to-noise ratio 'standard' profile. The standard profile 
is created by summing together all previous observa- 
tio ns of the pulsar using a t e chniqu e described in detail 
in IWeltevrede fc JohnstonI l|2008bl ). The topocentric 
TOA at the weighted centre frequency of the observa- 
tion is then added to the database of previous TOAs 
obtained for this pulsar. 

Pulsar timing is car ried out using the tempo2 pack- 
age (|Hobbs et al.l20(5^ ). In addition to other auxiliary 
files, TEMP02 requires a model for the spin behaviour 
of the pulsar and the database of TOAs. After con- 
verting the TOAs to t he solar system barycen tre us- 
ing t he DE405 model ij Edwards et al.1 2006: S tandishI 

11998'), TEMP02 compares the TOAs to the model pre- 
dictions to produce a set of residuals. These residuals 
are then used to refine the initial input model. 

Assuming we have an accurate pulsar position ob- 
tained through other means, the first step to deter- 
mining a timing solution is to fit only for the pulsar 
spin frequency, v, and frequency derivative, i>. This ef- 
fectively removes linear and quadratic terms from the 
residuals; what remains is then generally dominated 
by a cubic term. In these young pulsars, this cubic 
term is generally many orders of magnitude larger than 
the frequency second derivative expected from a simple 
dipole braking model and is the result of timing noise 
or glitch recovery (e.g. Johnston & Galloway 1999). 
Although one can, in principle, continue to fit higher 

-ocder polynomials to the data in order to whiten them, 

i-#iaifferent technique is used. 
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Table 1: Parameters and post-fit rms in milli-periods (mP) for the pulsars taken from Smith et al. (2008) 
observed at Parkes in support of the Fermi mission. 



Name 


P 


DM 


oL J 


rms 


Name 




P 


DM 


o L J 


rms 




fms) 


fpc/cm^) 


(erg/s) 

o/ ^/ 


(mP) 






fms) 


fpc/cm'') 


(crg/s) 


(mP) 


J0543+2329 


246.0 


78 


34.6 


1.1 


J1302- 


-6350 


47.8 


147 


35.9 


2.8 


J0614+2229 


335.0 


97 


34.8 


3.3 


J1305- 


-6203 


427.8 


470 


34.2 


2.0 


J0627+0705 


475.9 


138 


34.0 


0.4 


J1320 


5359 


279.7 


98 


34.2 


0.5 


J0659+1414 


384.9 


14 


34.6 


3.2 


J1327- 


-6400 


280.7 


681 


34.7 


8.1 


J0729-1448 


251.7 


92 


35.4 


5.3 


J1341- 


-6220 


193.3 


717 


36.1 


3.7 


J0742-2822 


166.8 


74 


35.1 


0.9 


J1349- 


-6130 


259.4 


285 


34.1 


0.9 


J0745-5353 


214.8 


122 


34.0 


4.2 


J1357- 


-6429 


166.1 


128 


36.5 


7.0 


J0821 3824 


124.8 


196 


34.7 


1.9 


J1359 


6038 


127.5 


294 


35.1 


0.2 


J0834-4159 


121.1 


240 


35.0 


3.9 


J 1406- 


-6121 


213.1 


542 


35.3 


9.4 


J0835-4510 


89.4 


68 


36.8 


2.3 


J1412- 


-6145 


315.2 


515 


35.1 


3.6 


J0855-4644 


64.7 


238 


36.0 


12.0 


J1413 


6141 


285.6 


677 


35.7 


89.7 


J0857-4424 


326.8 


184 


34.4 


1.4 


J1420- 


-6048 


68.2 


360 


37.0 


16.7 


J090 1-4624 


442.0 


199 


34.6 


0.6 


J1452- 


-5851 


386.6 


262 


34.5 


2.2 


J0905-5127 


346.3 


196 


34.4 


0.4 


J1452 


6036 


155.0 


350 


34.2 


0.4 


J0908 4913 


106.8 


180 


35.7 


0.5 


J1453 


6413 


179.5 


71 


34.3 


0.6 


J0940 5428 


87.5 


134 


36.3 


4.5 


J1509- 


-5850 


88.9 


138 


35.7 


12.1 


J0954-5430 


472.8 


200 


34.2 


0.7 


J1512- 


-5759 


128.7 


629 


35.1 


0.9 


J 1003-4747 


307.1 


98 


34.5 


0.5 


J1513- 


-5908 


150.7 


252 


37.3 


11.2 


J1015 5719 


139.9 


279 


35.9 


4.8 


J1514- 


-5925 


148.8 


194 


34.5 


3.2 


J1016-5819 


87.8 


252 


34.6 


2.0 


J1515- 


-5720 


286.6 


482 


34.0 


1.5 


J1016-5857 


107.4 


394 


36.4 


3.2 


J1524- 


-5625 


78.2 


153 


36.5 


5.3 


J1019-5749 


162.5 


1039 


35.3 


8.0 


J1524 


-5706 


1116.0 


833 


34.0 


0.6 


J1020 6026 


140.5 


445 


35.0 


12.9 


J 1530 


5327 


279.0 


50 


33.9 


0.8 


J1043 6116 


288.6 


449 


34.2 


0.5 


J1531 


5610 


84.2 


111 


36.0 


1.5 


J1048-5832 


123.7 


129 


36.3 


5.3 


J1538- 


-5551 


104.7 


603 


35.0 


8.4 


J1052-5954 


180.6 


491 


35.1 


12.1 


J1539- 


-5626 


243.4 


176 


34.1 


0.9 


J1055 6032 


99.7 


633 


36.0 




J1541- 


-5535 


295.8 


428 


35.0 


5.3 


J1057-5226 


197.1 


30 


34.5 


0.3 


J1543- 


-5459 


377.1 


346 


34.6 


1.9 


J1105-6107 


63.2 


271 


36.4 


5.2 


J1548- 


-5607 


170.9 


316 


34.9 


1.8 


J1112 6103 


65.0 


599 


36.7 


10.5 


J1549 


4848 


288.3 


56 


34.4 


1.2 


J1115-6052 


259.8 


228 


34.2 


0.7 


J1551- 


-5310 


453.4 


493 


34.9 


7.1 


J1119-6127 


407.7 


707 


36.4 


0.2 


J1600- 


-5044 


192.6 


261 


34.4 


0.4 


J1123-6259 


271.4 


223 


34.0 


0.9 


J1600- 


-5751 


194.5 


177 


34.0 


0.9 


J1124 5916 


135.3 


330 


37.1 




J1601 


5335 


288.5 


195 


35.0 


9.4 


J1138-6207 


117.6 


520 


35.5 


4.9 


J1611- 


-5209 


182.5 


128 


34.5 


0.3 


J1156-5707 


288.4 


244 


34.6 


0.7 


J1614- 


-5048 


231.7 


583 


36.2 


5.6 


J1216-6223 


374.0 


787 


34.1 


4.4 


J1617- 


-5055 


69.4 


467 


37.2 




J1224 6407 


216.5 


97 


34.3 


0.3 


J1626- 


-4807 


293.9 


817 


34.4 


17.8 


J1248 6344 


198.3 


433 


34,9 


8.0 


.T1627 


4706 


140.7 


456 


34,4 


9.4 


J1301 6305 


184..5 


374 


36.2 


14.0 


.J1632 


4757 


228. G 


578 


34,7 


7.8 
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Table 1: Parameters and post-fit rms in milli-pcriods (mP) for the pulsars taken from Smith et al. (2008) 
observed at Parkes in support of the Fermi mission (cont). 



Name 




P 


DM 


log\E] 

O L J 


rms 


Name 


P 


DM 


O L J 


rms 






fms) 


(pc/cm^) 


(erg/s) 

o/ 


(mP) 




(ms) 


(pc/cm^) 


(erg/s) 

o/ 


(mP) 


J1632- 


-4818 


813.5 


758 


34.7 


12.8 


J 180 1-2451 


124.9 


289 


36.4 


4.9 


J1637 


4553 


118.8 


193 


34.9 


0.9 


J1803-2137 


133.6 


234 


36.3 


6.5 


J1637~ 


-4642 


154.0 


417 


35.8 


9.9 


J1806-2125 


481.8 


750 


34.6 


3.8 


J 1638- 


-4417 


117.8 


436 


34.6 


4.0 


J1809-1917 


82.7 


197 


36.3 


5.3 


J1638- 


-4608 


278.1 


424 


35.0 


0.9 


J1812-1910 


431.0 


892 


34.3 


9.0 


J1640- 


-4715 


517.4 


592 


34.1 


3.7 


J1815-1738 


198.4 


728 


35.6 


3.3 


J1643- 


-4505 


237.4 


484 


35.0 


2.4 


J1820 1529 


333.2 


772 


34.6 


5.5 


J1646- 


-4346 


231.6 


490 


35.6 


6.2 


J1824-1945 


189.3 


225 


34.5 


0.3 


J1648- 


-4611 


165.0 


393 


35.3 


6.2 


J1825-1446 


279.2 


357 


34.6 


1.7 


J1649 


4653 


557.0 


332 


34.0 


1.6 


J1826 1334 


101.5 


231 


36.4 


7.7 


J1650- 


-4502 


380.9 


320 


34.0 


2.0 


J1828-1057 


246.3 


245 


34.7 


6.7 


J1650- 


-4921 


156.4 


230 


34.3 


1.0 


J1828-1101 


72.1 


607 


36.2 


5.1 


J1702- 


-4128 


182.1 


367 


35.5 


6.5 


J1830-1059 


405.0 


162 


34.6 


1.0 


J1702 


4305 


215.5 


538 


34.6 


5.9 


J1831 0952 


67.3 


247 


36.0 


6.3 


J 1702 


-4310 


240.5 


377 


35.8 


4.1 


J1832-0827 


647.3 


301 


34.0 


0.7 


J 1705- 


-3950 


318.9 


207 


34.9 


2.0 


J 1833-082 7 


85.3 


411 


35.8 


1.7 


J1709- 


-4429 


102.5 


76 


36.5 


4.6 


J 1834-0731 


513.0 


295 


34.2 


4.3 


J1715- 


-3903 


278.5 


313 


34.8 


5.7 


J1835-0643 


305.8 


473 


34.7 


4.6 


J1718- 


-3825 


74.7 


247 


36.1 


3.0 


J1835-0944 


145.3 


277 


34.7 


5.9 


J1721- 


-3532 


280.4 


496 


34.7 


1.0 


J1835-1106 


165.9 


133 


35.3 


2.4 


J1722- 


-3712 


236.2 


100 


34.5 


1.7 


J1837-0559 


201.1 


318 


34.2 


4.4 


J1723 


3659 


202.7 


254 


34.6 


1.1 


J1837 0604 


96.3 


462 


36.3 


17.1 


J1726 


3530 


1110.3 


727 


34.5 


6.3 


J1838 0453 


380.8 


621 


34.9 


2.8 


J 1730- 


-3350 


139.5 


259 


36.1 


1.3 


J1838-0549 


235.3 


274 


35.0 


4.3 


J1731- 


-4744 


829.8 


123 


34.0 


0.4 


J1839-0321 


238.8 


449 


34.6 


4.1 


J1733 


3716 


337.6 


154 


34.2 


1.0 


J1839-0905 


419.0 


348 


34.1 


1.9 


J 1734- 


-3333 


1169.2 


578 


34.7 


9.6 


J1841-0425 


186.1 


325 


34.6 


0.5 


J1735- 


-3258 


351.0 


754 


34.4 


14.9 


J1841-0524 


445.7 


289 


35.0 


4.6 


J1737 


3137 


450.4 


488 


34.8 


4.2 


J1842 0905 


344.6 


343 


34.0 


1.1 


J 1738- 


-2955 


443.4 


223 


34.6 


0.9 


J1843-0355 


132.3 


798 


34.3 


12.4 


J1739- 


-2903 


322.9 


139 


34.0 


0.3 


J1843-0702 


191.6 


228 


34.1 


2.6 


J1739- 


-3023 


114.4 


170 


35.5 


5.8 


J1844 0256 


273.0 


820 


34.7 


23.7 


J1740 


3015 


606.8 


152 


34.9 


9.6 


J1844 0538 


255.7 


413 


34.4 


1.1 


J1745- 


-3040 


367.4 


88 


33.9 


0.7 


J1845-0743 


104.7 


281 


34.1 


0.3 


J1756- 


-2225 


405.0 


326 


34.5 


1.1 


J1847-0402 


597.8 


142 


34.0 


0.7 


J1757- 


-2421 


234.1 


179 


34.6 


0.4 


J1853+0011 


397.9 


569 


34.3 


2.9 


J1801- 


-2154 


375.3 


388 


34.1 


3.1 


J1853-0004 


101.4 


438 


35.3 


1.6 


J1801- 


-2304 


415.8 


1074 


34.8 


4.0 


J1903-h0925 


357.2 


162 


30.7 
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Table 2: Parameters and post-fit rms in milli- 
periods (mP) for tlie pulsars observed in addition 
to those contained in Smith et al. (2008). 



Table 3: Parameters of the 14 glitches observed in 
the timing data. 



^1 O TTl/i 

i \ cLllit: 




P 




loo-[Fl 
iOg[£/J 


rms 






(ms) 


VPC/cm ) 


lerg/sj 




Tfll HQ 
JUlUo" 


1 /I Q1 


oU ( .0 





oU.o 


u.o 


Tfi/ini 
JU4U1- 


- ( OUo 


040. 


00 


QO 
oZ.O 


1 n 




- 1 04o 


1 9/1 r; Q 


1 8 
lo 


ol.l 


n /I 

U.4: 


10630- 


-2834 


1244.4 


34 


32.2 


1.2 


J0738- 


-4042 


374.9 


161 


33.1 


0.6 


J1028- 


-5819 


91.4 


97 


35.9 


0.4 


J1456- 


-6843 


263.4 


9 


32.3 


15.9 


J1602- 


-5100 


864.2 


171 


33.6 


0.9 


J1638- 


-4725 


763.9 


550 


32.4 


52.1 


J1705- 


-1906 


299 


23 


33.8 


0.3 


J1825- 


-0935 


769 


19 


33.7 


0.9 


J1845- 


-0434 


486.8 


231 


33.6 


0.8 



The technique (developed by Hobbs et al. 2004) in- 
stead removes harmonically related sinusoids from the 
residual data after fitting for v and v and continues to 
remove components until the residuals resemble white 
noise. This whitening technique is implemented using 
the FITWAVES algorithm within tempo2. 

An added complication arises if the pulsar suffers 
a glitch. A glitch causes an abrupt change in both v 
and ii and the magnitude of the change can be large. 
With sparsely sampled data it is very difficult to main- 
tain phase coherence through the epoch of the glitch. 
This implies that the glitch epoch cannot be well de- 
termined - in the case of our data, we generally cannot 
determine the glitch epoch to better than ~15 days. 
In order to measure the glitch parameters, we deter- 
mine a timing solution before and after the glitch and 
then, using the two solutions, estimate Ai/ and A!>. 
With these estimates, and by setting the glitch epoch 
midway between the observation just before and just 
after the glitch, a timing solution can be obtained for 
the entire data set with the addition of a phase jump 
at the glitch epoch. 

One of the important parameters to measure, in 
the context of models of gamma-ray emission, is the 
offset between the radio and gamma-ray profiles. In or- 
der to do this, an accurate value of the dispersion mea- 
sure (DM) is needed so that the dispersion delay be- 
tween the ~1 GHz radio profile and the ~GeV gamma- 
ray profile can be corrected for. We measured the DM 
across the 256 MHz bandwidth at 1.4 GHz assuming 
no profile evolution across the band which is generally 
sufficient to determine the DM to ~0.01 cm~^pc. This 
corresponds to an uncertainty of ~40 ^s in the offset 
between the radio and gamma-ray profiles. 



4 Results - General Timing 

In the last columns of Tables [1] and [2] we list the rms of 
the residuals in milli-periods for each pulsar following 



Name 




Epoch 










(MJD) 






J0729- 


-14481' 


54711('21i 


6 6203f4i 


13 86f5i 


J1048- 


-5832 


54495(10) 


3.0431(2) 


4 77f2i 


J1052- 


-5954'l' 


54526(22) 


6.7478(7) 


14 1(A) 


J1105- 


-6107 


54711(21) 


0.023(1) 


0.0(3) 


J1341- 


-6220 


54468(18) 


3.0782(3) 








54870(11) 


3.0667(3) 




J1357- 


-6429 


54803(17) 


1.752(7) 


2.8(2.8) 


J1410- 


-6132+ 


54652(19) 


0.2726(5) 


-3.31(5) 


J1420- 


-6048t 


54652(20) 


0.9371(3) 


2.95(1) 


J1709- 


-4429 


54710(22) 


2.7497(1) 


4.95(1) 


J1737- 


-3137t 


54352.3(1.0) 


1.3425(4) 


0.4(6) 


J1740- 


-3015t 


54468(17) 


0.037(1) 


-0.4(4) 


J1801- 


-2451 


54680(9) 


3.0904(2) 


5.04(2) 


J1841- 


-0524t 


54495(10) 


1.0359(4) 


0.11(15) 



TPulsar not previously known to have glitched. 



fitting for u, v and whitening if necessary. When creat- 
ing gamma-ray profiles from the radio ephemeris, the 
rms constrains the maximum useful number of phase 
bins across the light curve. The rms is less than 2 milli- 
periods for 64 pulsars in our sample - 500 bins across 
the light curve is then possible. Only for 18 pulsars is 
the rms worse than 10 milli-periods. We note that the 
accuracy of the time tagging of the gamma-ray photons 
is significantly better than 1 /is (Smith et al. 2008) 
and is not a limiting factor in the time-resolution of 
the resulting gamma-ray profiles. Typically, the small 
number of photons received from gamma-ray pulsars 
limits the profile resolution to less than 100 bins. 

In three cases we have been unable to maintain 
phase coherence in the timing data. The pulsars con- 
cerned are PSRs J1055-6032, J1617-5055 and J1903-h0925. 
In the latter case we have less than 1 yr of timing data 
but the first two named pulsars appear to have ex- 
tremely large short-term timing noise. They need to 
be monitored on a much more regular basis than once 
per month in order to maintain phase coherence. In 
additio n, PSR J1124-5916 is an extremely weak radio 
pulsar (jCamilo et al.l [200^ 1 for which it takes ~5 hr to 
obtain a timing point using Parkes and it is therefore 
not observed as part of our timing program me. Follow- 
ing its detection in gamma-rays by Fermi (jAbdo et al.l 
l2009d '). we observed it on a single occasion in order to 
phase align the radio and gamma-ray profiles. 



5 Results - Glitches 

We have detected 14 glitches in 13 of the pulsars in our 
sample. These glitches are tabulated in Table [3] with 
the parameters estimated using the scheme outlined in 
Section 3. The epoch of the glitch is set at the mid- 
point between the last pre-glitch and first post-glitch 
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Figure 1: Timing residuals for PSRs J0659+1414 (top left), J0742~2822 (top right), J0835-4510 (bottom 
left) and J1028-5819 (bottom right). Residuals before (top) and after (bottom) the use of fitwaves are 
shown for the latter three pulsars. Note the very different scales on the y-axis for these cases. 



data point. We note that Fermi data can be used to 
better c onstrain the glitch epoch for PSR J1709-4429 
(|Abdo e t al. 2009c). The third column of the table 
gives the fractional change in the spin frequency, Av/i/ 
and the final column lists the fractional change in the 
frequency derivative, Ai>/i>. The listed errors are twice 
the formal errors given by tempo2 and given in brack- 
ets which refer to the last digit(s). Three of the pulsars 
in the table have been detected by Fermi. For PSRs 
J1048-5832 and J1420-6048 the glitch occurred prior 
to the launch of the satellite. The glitch in PSR J 1709- 
4429 happened shortly after Fermi began its all sky 
survey in 2008 August. 

Of the 13 pulsars, 7 have not previously been known 
to glitch, these are listed with the f symbol in Ta- 
ble 13] The other 6 pulsars have previously been ob- 
served to g l itch a s tabulated in the compilation of 
IWang et al.l (|2000l '). Most notable is PSR J1341-6220 
(B1338-62) which has the shortest averag e interval be- 
tween glitches of any pulsar, estimated bv lWang et al.l 
(I2OOOI ') to be 250 days. In the 800 days of our tim- 
ing data we see two glitches separated by some 400 
days. Both glitches have about the same magnitude 
and are larger than any of the other 12 glitches re- 
ported for this pulsar. There is not sufficient time 
between the glitches to determine a reliable value of 



A!>/i>. PSR J1740-301 5 (B1737-30) is another pul- 
sar which glitches often ijZou et al ] |2008 ll and we have 
detected one relatively small glitch in this object. 

The smallest glitch we report here has Av/u of 2 x 
10~* and we believe our sample to be complete above 
a limit of about 10~*. Below this limit it is very hard, 
with our data set, to distinguish between glitches and 
timing noise. We attempted to fit glitch parameters 
to possi ble 'cuspy' looking r e sidua ls, similar to those 
seen in Ijanssen fc Sta ppersi (j200d ) for example, but 
we find that fitting glitches and using fitwaves return 
very similar results in spite of their different functional 
forms. 

6 Results - Fermi pulsars 

In this section we describe in detail the ten pulsars 
which have been confirmed as gamma-ray emitters by 
Fermi (|Abdo et al ] l2009d ). We show residuals as a 
function of time for each pulsar and provide the com- 
plete ephemerides in Tables |4] and [5] including details 
of the FITWAVES parameters where appropriate. These 
ephemerides, in conjunction with the publicly avail- 
able photon list from Fermi, can be used to generate 
gamma-ray profiles of the Ferroi-detected pulsars. 



www. publish . csiro. au /journals /pasa 



7 



PSR .J1048-5B33 


• • • 


: 

, 1 


H \ ' \ ' h 

• • 
■ ^ 

• • • 

* 



PSR Ji057-5336 



li 




600 800 
since MJD 54000 



PSR Ji509-5650 



4 . + 



400 600 

Days since MJD 54000 



600 800 
Days since XIJD 54000 



Figure 2: Timing residuals for PSRs J1048-5832 (top left), J1057-5226 (top right), J1420-6048 (bottom 
left) and J1509-5850 (bottom right). Residuals before (top) and after (bottom) the use of fitwaves are 
shown for the latter three pulsars. Note the very different scales on the y-axis for these cases. The arrow 
indicates the epoch of the ghtch for PSRs J1048-5832 and J1420-6048. 



We note that for all the pulsars, the clock correc- 
tion is to TT(TAI), the solar system ephemeris used 
is DE405, and we are using the barycentric dynamical 
time (BDT) definition for the spin parameters. For a 
full description of the meaning of these terms and their 
implementation in tempo2 see Hobbs et al. (2006). 

The pulsar's position, proper motion and parallax 
are held constant during the fitting procedure using 
the best values taken from the literature. The pulsar's 
DM is measured across the band at 1.4 GHz and is 
also held constant during the fitting. The epoch refers 
to the epoch at which the rotational frequency and po- 
sitional parameters are measured. The reference MJD 
and reference frequency given in the table is the TOA 
of the fiducial point (the peak of the pulse profile) at 
the Parkes telescope. Start and finish MJD denote the 
first and last data points used and hence the range of 
dates over which the fit is valid. The accuracy of the 
ephemerides outside the valid MJD range deteriorates 
quickly in the majority of cases. 

PSR J0659+1414/B0656+14 (seeFig[l]): This 
pulsar is a nearby, relatively low E pulsar. It is known 
to pr oduce both thermal and non-thermal X-ray emis- 
sion (|d7 Luca et al.ll2005l ) and suspected as a gamma- 



ray pulsator (jRamanamurthv et al.|[T99 6*) before con- 
firmation with Fermi (Weltevrede et al. 2009). Bright, 
"spiky" emission is occas ionally seen in the radio pulses 
(|Weltevrede et al]|2006t ). The pulsar has low timing 
noise and is not known to have glitched in spite of 
decades of monitoring. There was no need to whiten 
the data using fitwaves as the timin g noise is low . 
Astrometric data are taken from lBrisken et al.l lj2003t ). 

PSR J0742-2822/B0740-28 (see Fig[T|): This 
pulsar was not known as a high-energy emitter un- 
til its recent detection with Fermi (Weltevrede et al. 
2009). It is extremely bright in the radio and the 
rotation axis a ppears to point in the direction of its 
proper motion (| Johnston et al.ll2005t) . The pulsar is 
not known to glitch but suffers a high level of timing 
noise which was successfully removed with the fit- 
WAVES algorithm. A strometric data are taken from 
iFomalont et all ([1993). 

PSR J0835-4510/B0833-45 (see Fig[l|): The 

Vela pulsar was one of the first pulsars discovered and 
has a timing history stretching back over 40 years. 
Located in the Vela supernova rem nant, its gamma- 
ray p ulsations have long been known (jThompson et ahl 
119751 ') with the EGRET results discussed in Kanbach 
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Figure 3: Timing residuals for PSRs J1709-4429 
(top), J 1718-3825 (bottom). Residuals before 
(top) and after (bottom) the use of fitwaves are 
shown. Note the very different scales on the y- 
axis. The arrow indicates the epoch of the glitch 
for PSR J1709-4429. 



et al. (1994) and recent results from the Agile satel- 
lite are reported in Pellizoni et al. (2009). Fermi has 

S iroduced beautiful high quality data on this pulsar 
Abdo et al] 12009^1. Vela suffer s glitches on a quasi- 
regular basis l|Wang et al1l2000l ). but no large glitch 
has occurred over the timespan of our data. The po- 
sitio n, proper moti o n and parallax for Vela are taken 
from lDodson et all ||2003D . 

PSR J1028-5819 (see Fig[l]): This pulsar was 
discovered in 2008 March in a search of unidentified 
EGRET error boxes ijKeith et al.ll2008l ) and subse quently 
detect ed as a pulsed gamma-ray emitter by Fermi (|Abdo et 
12009^ 1. The very narrow radio profile allows for ac- 
curate TOAs to be obtained and the residual is only 
40/is. The puls ar position is from the interferometric 
observations of iKeith etHl (|2008h . 

PSR J1048-5832/B1046-58 (seeFig[2D: This 
pulsar is a 'Vela-like' object, proposed tentatively as a 
gamm a-ray emitter from EGRET data (Ka spi et al.l 
Hood ) and confirmed by Fermi (jAbdo et all l2009bD . 
Large pulses are occasionally seen on th e rising edge 
of the profile (j Johnston fc Romanill200j ). The pulsar 
glitched in early 2008. The timing noise is large and 



tracks smoothly ove r the glitch e poch. The pulsar's 
position is that from lWang et all (|2000h . 

PSR J1057-5226/B1055-52 (seeFiglH): This 
pulsar has an interpulse and both the main and inter- 
pulses have a co m plex p ulse morphology, described ex- 
tensively in I Weltevrede fc Wright (2009). Although it 
has a relatively low E, it was seen as a pulsed gamma- 
ray emitter by EGRET (Fierro et al. 1993). The 
gamma-ray emission is claimed to be associated with 
the pole that produces the radio interpulse emission 
(e.g. Weltevrede & Wright 2009). The timing noise 
seems quasi-periodic with a period close to one year. 
The pulsar's position is taken from the radio timing 
paper of Newton et al. (1981). 

PSR J1420-6048 (see Fig [2]): This pulsar lies 
in a complex region of th e Galactic plane and is a 
known X-ray emitter ( Robe rts et al ]'2001). It was 
mooted as a p ossible counterpart to an EGRET uniden- 
tified source (|D'Amico fc et al.ll200lh . and confirmed 
as a gamma-ray emitter by Fermi (Weltevrede et al. 
2009). A glitch in the pulsar happened in mid 2008, 
just before Fermi began its all sky survey. The timing 
noise is large. The pulsar's position is from D'Amico 
et al. (2001). 

PSR J1509-5850 (see Fig [g]): The pul sar is 
locate d inside an X-rav ijKargaltsev et aLll2008D and 
radio (|Hui fc Becked 120071 )" wind nebula. Weltevrede 
et al. (2009) report the gamma-ray detection of the 
pulsar. Rs timing noise is very low, the residual after 
fitting for u and i/ is only 20 milli-periods and there is 
no need for whitening using fitwaves. The pulsar's 
position is taken from the discovery paper (Kramer et 
al. 2003). 

PSR J1709-4429/B1706-44 (seeFiglS]): This 
pulsar is a 'Vela-like' object, detected as a gamma-ray 
emitter by EGRET (Thompson et al. 1992). A glitch 
occurred in 2008 August; the daily observations and 
monitorin g by Fermi allows an accurate epoch to be 
obtained (I Abdo et al.ll20093 ). The timing noise contin- 
ues smoothl y through the glitc h. The pulsar position 
is that from lWang et all l|2000"'). 

PSR J1718-3825 (see Fig[3l): The pulsa r has 
an associated X-ray nebula llHinton et al.l 120071 ) and 
a TeV source (Aharonian et al. 2007). The gamma- 
ray detection of the pulsar is reported in Weltevrede 
et al. (2009). Its timing noise is relatively large and 
has a longer timescale than many of the other pulsars 
in our sample. The pulsar's position is taken from 
Manchester et al. (2001). 

Summary 

We started a comprehensive program of pulsar timing 
with the Parkes radio telescope in 2007 February. The 
launch of the Fermi satellite in 2008 June has resulted 
in the detection of more than 20 (non-MSP) radio pul- 
sars at gamma-ray wavelengths, of which 10 are ob- 
served at Parkes. In this paper we have provided the 
timing ephemerides for the 10 pulsars, which, in con- 
junction with the publicly available photon parameters 
from Fermi, allow a recreation of the gamma-ray pro- 
files. In addition we have detected 14 glitches in 13 
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difFerent pulsars and provide the glitch parameters. 

The Fermi mission will continue for at least a fur- 
ther 4 years and we plan to carry out the pulsar timing 
for a similar length of time. Researchers who require 
ephemerides from any of the pulsars contained in Ta- 
ble [T] can contact the authors of this paper. 
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Table 4: Timing parameters for 5 pulsars detected in gamma-rays by Fermi. 



Param 



J0659+1414 



J0742-2822 



J0835-4510 



J1028-5819 



J1048-5832 



Right Ascension (J2000) 
Declination (J2000) 

Spin Frequency, u, {s~^) 
y (10-15) 
DM (cm-3pc) 

Proper Motion in RA (mas/yr) 

Proper Motion in Dec (mas/yr) 

Epoch (MJD) 

Ref. MJD 

Ref. Frcq (MHz) 

Glitch Epoch (MJD) 

Glitch Phase Offset 

Glitch Ai^ (10-5) 

Glitch (10-1^) 

Start (MJD) 

Finish (MJD) 

Wave Frequency (r ad/day) 

Wave 1 

Wave 2 

Wave 3 

Wave 4 

Wave 5 

Wave 6 

Wave 7 



06:59:48.134 
+14:14:21.5 

2.598136794(2) 
-370.981(5) 
13.7(2) 
44.07 
-2.4 
49721.0 
54762.6946809825 
1369.235 



54505 
55011 



07:42:49.026 
-28:22:43.70 
5.996283159(8) 
-611(7) 
73.790(3) 
-29 
4 

54615.0 
54632.9827544637 
1377.652 



54220 
55010 
0.0063560 
-0.4498,+0.05777 
+0.09547, 0.02979 
-0.03240,+0.0167 
+0.01348,-0.009575 
0.004789,+0.005083 
+0.001260,-0.002011 
-0.0003854,+0.0008335 



08:35:20.61149 
-45:10:34.8751 
11.191119(5) 
-15600(170) 
67.9754(8) 
-49.68 

29.9 
54217.0 
54606.2500568554 
1373.864 



54207 
55010 
0.0052116 
+3.318. 1.538 
+0.18G3, 1.021 
-0.3375,-0.2816 
-0.1746,+0.08070 
+0.003159,+0.07619 
+0.02241,+0.007917 
+0.003275,-0.003641 



10:28:27.95 
-58:19:05.225 
10.9404949958(3) 
-1928.4(3) 
96.39(9) 



54787.0 

54786.6818285185 
1374.656 



54563 
55011 
0.0093621 
0.008672, 0.002189 
+0.002582, 0.001674 
-0.0002706,-3.2135e-05 



10:48:12.2 
-58:32:05.8 
8.08398924(6) 
-6259(7) 
128.822(8) 



54615.0 
54606.2571751136 
1373.864 

54490 
-0.484(6) 
2.4616(3) 
-3.7(4) 
54220 
55011 
0.0059587 
0.3531, 0.1548 
+0.07018,+0.01160 
-0.01741, 0.004765 
+0.007688,+0.0005985 
-0.001823,+0.002077 



Table 5: Timing parameters for 5 pulsars detected in gamma-rays by Fermi. 



Param 


J1057-5226 


J1420-6048 


J1509-5850 


J1709-4429 


J1718-3825 


Right Ascension (J2000) 


10:57:58.84 


14:20:08.237 


15:09:27.13 


17:09:42.728 


17:18:13.565 


Declination (J2000) 


-52:26:56.3 


-60:48:16.43 


-58:50:56.1 


-44:29:08.24 


-38:25:18.06 


Spin Frequency, v, (s~^) 


5.0732274494(7) 


14.66248295(6) 


11.2455254698(1) 


9.756518994(9) 


13.391571381(4) 


!> (lO-l-"^) 


-150.4(6) 


-17710(16) 


-1159.23(1) 


-8831.0(7) 


-2375(4) 


Epoch (MJD) 


54615.0 


54615.0 


54615.0 


54616.0 


54616.0 


DM (cm-3pc) 


30.1 


360 


137.7 


75.69 


247.4 


Ref. MJD 


54606.2654505538 


54633.2224514535 


54633.2560188871 


54633.4457155762 


54633.4498848905 


Ref. Preq (MHz) 


1373.864 


1377.652 


1377.652 


1377.652 


1377.652 


Glitch Epoch (MJD) 






54660 




54692.8 


Glitch Phase Offset 






-0.23(2) 




-0.178(9) 


Glitch Aiy (10-5) 






1.3659(4) 




2.6837(3) 


Glitch Az> (10-") 






-10.0(3) 




6.2(1) 


Start (MJD) 


54220 


54220 


54220 


54220 


54220 


Finish (MJD) 


55011 


55011 


55011 


55011 


55011 


Wave Frequency (rad/day) 


0.0059587 


0.0039726 




0.0061797 


0.0059593 


Wave 1 


-0.02318,-0.0007458 


+5.079,-1.557 




+0.06933,-0.2996 


-0.3333,+0.007675 


Wave 2 


+0.005196,+0.0001067 


-0.6315, +0.04927 




+0.01173,+0.0304 


+0.06940,-0.009321 


Wave 3 


+0.0009044 ,+0.0004835 


+0.08892,-0.007211 




-0.002489,-0.001332 


-0.02062,+0.006114 


Wave 4 


+0.0003618,+0.0003014 








+0.005628,-0.003100 


Wave 5 


-0.0001134,-2.6808e-05 








-0.0009510,+0.0008923 



